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ДЕЯКІ АНАЛІТИЧНІ АПРОКСИМАЦІЇ 
ДЛЯ НАБЛИЖЕНОГО ВИЗНАЧЕННЯ КОНТАКТУ 


Вказано спосіб заміни нелінійної контактної задачі гравіметрії на лінійну. Подано аналітичний ви- 
раз для її розв'язання з прийнятною точністю. Виявлено, що відомі аналітичні конструкції для обчислен- 
ня густинних контактів є його частинними випадками. Наводяться ряд ітераційних поправок для уточне- 
ння наближень, які є аналітичними апроксимаціями шуканого контакту. Вказана міра уточнення цих наб- 


лижень. 
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Вступ. У процесі становлення нової парадигми 
в теорії інтерпретації потенціальних полів |Стра- 
хов, 2001) актуальною залишається розробка но- 
вого математичного апарату для моделювання по- 
тенціальних геофізичних полів, який був би адек- 
ватним сучасній геофізичній практиці. Про адек- 
ватні вимоги до нового математичного діалекту 
геофізичного моделювання та ключову роль аналі- 
тичних апроксимацій геологічного середовища і 
геофізичного поля сказано в роботах Страхова, зо- 
крема |Страхов, 2007). Продовжуючи наявні роз- 
робки у цій сфері, пропонуємо нові аналітичні ап- 
роксимації для горизонтально-шаруватого геоло- 
гічного середовища, що мають місце в рамках під- 
класу відомого класу контактних поверхонь Стра- 
хова. 


До постановки проблеми. Розв'язок сх) не- 
лінійної оберненої задачі потенціалу для контакт- 
ної поверхні |Нумеров, 1930| на класі Нумерова 

1, и - 
Ми! (ап) меж сх) поділу однорідних сере- 
довищ з високою точністю можна замінити на лі- 


нійне рівняння, яке характеризує аналітичне про- 
довження вертикальної складової напруженості 


поля и(х) до тяжіючих мас на середню глибину 
па сх) | Старостенко и др., 19881: 
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Аналізуючи фізичну суть процедури (1) при обме- 








й І, , 
женнях, які властиві саме класу Ми! 1, п), пі- 


сля ряду перетворень отримуємо таке лінеаризова- 
не наближення 


ІЗ 
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Вираз (2) апроксимує розв'язок задачі з точністю 
. й коїч ж -р 
до квадрата різниці значень поля Й -й б 


на рівнях у -« 0 й у - -/. Проте цей вираз не 


дає змоги визначити поле их, п) , просто продов- 
живши и(х) на рівень у - й » 0. Адже пряма 
уз Й розтинає тяжіючі маси надвоє і одна з час- 
тин при с(х)х Й буде нижча рівня у - П. Нас- 
лідком цього геометричного факту є те, що функ- 
ція их, у) при у « Й - негармонічна. 
Розв'язання задачі. Однак знайти наближення 


сх, 1) можна не прямо за виразом (2), а опосере- 


дковано - як межу послідовності с, (п), роз- 


в'язків, які отримують за допомогою процесу іте- 
рацій типу Андрєєва 
17 й 
Савь Нє б (5 п)- Ї 2 2 Ех 
пй (6-х) 


м(х), сп) аз м) і и(х)- паб, (3) 
з наступною регуляризацією за принципом Лавре- 


нтьєва (цього цілком достатньо для отримання пе- 
ршого наближення). Доведено, що на класі 





1 зла щі 
Ми! 4) (п) розв'язок нелінійної контактної за- 


дачі сх) із точністю до в(во-а) при 


Й -500 обчислюється за допомогою такого лінеа- 


ризованого виразу 
1 пі 
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З цієї формули, як окремі випадки при 
п - 0,1,2 , випливають відомі аналітичні констру- 


кції |Нумеров, 1930), |Ниєрез, 1942| 1 |Страхов, 


19761, де Аи(2) - ц(Є,0)-- и(х, 0): 
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Цікаво, що формули (4) та (б), (7) при п2 3 
неефективні для обчислення наближень контакту 
- через наростання з номером ітерації похибок ок- 
руглень, особливо при обчисленні в околі точки 





5 з х.. Виявлено, що цей ефект прогресії похибок 
не властивий для виразу (5). 

Щоб уникнути такого небажаного накопичення 
похибок, вдаємось до обхідного маневру - ядра ін- 


іо (5-Й 


тегралів підсумовуємо з різними знаками у напря- 
мку зростання біномних коефіцієнтів. Практика 
обчислень на простих чисельних моделях виявила, 
що підсумування елементарних дробів у ядрах пе- 
ретворень з індексами п » 0 за схемою (4) визна- 
чає додатні ядра перетворень лише при п - 1,2,3, 


апри по4 - уже від'ємні, що свідчить про роз- 

біжність відповідних інтегральних квадратур. 
Тому за результатами чисельного моделюван- 

ня виявлено, що для практичних обчислень набли- 


жень (4) контакту 0/ (п) (п) достатньо узяти за 


"нульове наближення" вираз (5) і обмежитись од- 
нією з формул: 
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ЗОВНІ 


Вибір детальності обчислення наближення за- 
лежить від характеру вхідного поля и(х) і масш- 
табу побудов (міри задачі). Слід зазначити справе- 
дливість вказаних вище конструкцій лише для об- 
ластей малої міри (4 100 км') і нескладної геомет- 
рії. 

Уточнення моделі контакту. Наближення (4) 
«РЯД; можна істотно уточнити за допомогою 
ряду поправок, які є аналітичними апроксимація- 
ми шуканого контакту різної міри детальності на 
зазначеному класі Ми) (1 п) шуканих розв'яз- 
ків. 

Наприклад, поправка, яка має конструкцію 
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- - -3 3 
підвищує точність наближень до иа (2). Одна- 





че, чисельне моделювання виявило, що в умовах 
слабоградієнтного поля достатньо дати поправку 
спрощеного вигляду 


АС, 1) ко (в, 


де (2) - деяке "середнє" (осереднене засобами 


удць-ю), п0) 


фільтрації) значення поля на поверхні спостере- 
жень. 
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льного рівняння 1-го й 
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Ця поправка (11) допускає узагальнення у виг- 


ляді чисельної конструкції б Х ЗХ з с Сн д6 (»), 
де Дб, (х) - йт Дт хх) отримують як ряд пос- 
тою 
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лідовних наближень за допомогою ітераційного 
процесу: 
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Висновки. Відтак, отримано ряд аналітичних 
конструкцій (8-12), застосування яких щодо відо- 
мих раніше апроксимаційних зображень контакту 
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Доведено, що на класі Мио(1, П) послідовні 


наближення (8-9) 1 (11-12) є однозначними і стій- 
кими (збіжними). Чисельно кожен спосіб зводить- 
ся до розв'язання лінійного інтегрального рівнян- 
ня Фредгольма 1-го роду з ядром типу Пуассона. 
У перспективі очікується аналіз точності запропо- 
нованих конструкцій відносно відомих методів 
|Гравиразведка.., 1981) та апробація на польових 
даних. 
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НЕКОТОРЬКЕ АНАЛИТИЧЕСКИК АППРОКСИМАЦИИ ДЛЯ 
ПРИБЛИЖЕННОГО ОПРЕДЕЛЕНИЯ КОНТАКТА 


Ю.Х-ЙИ. Дубовенко 


Указан способ заменьт нелинейной контактной задачи гравиметрии линейной. Дано аналитическоє 
вьтражение для ее решения с приемлемой точностью. Вьіявлено, что известньєе аналитические конструк- 
ции для вьтчисления плотностньх контактов есть его частньми случаями. Наводится ряд итерационньїх 
поправок для уточнения приближений, каковьте являются аналитическими аппроксимациями искомого 
контакта. Указана мера уточнения зтих приближений. 
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